Decoherence-Free Subspaces for Multiple-Qubit Errors: (I) Characterization 

Daniel A. Lidar, 1 Dave Bacon, 1:2 Julia Kempe 1,3 ' 4 and K.B. Whaley 1 
Departments of Chemistry 1 , Physics 2 and Mathematics 3 , 
University of California, Berkeley, CA 94720 
Ecole Nationale Superieure des Telecommunications, 
Paris, France 
(February 1, 2008) 

Coherence in an open quantum system is degraded through its interaction with a bath. This 
decoherence can be avoided by restricting the dynamics of the system to special decoherence-free 
subspaces. These subspaces are usually constructed under the assumption of spatially symmetric 
system-bath coupling. Here we show that decoherence-free subspaces may appear without spatial 
symmetry. Instead, we consider a model of system-bath interactions in which to first order only 
multiple-qubit coupling to the bath is present, with single-qubit system-bath coupling absent. We 
derive necessary and sufficient conditions for the appearance of decoherence-free states in this model, 
and give a number of examples. In a sequel paper we show how to perform universal and fault 
tolerant quantum computation on the decoherence-free subspaces considered in this paper. 
PACS numbers: 03.67.Lx, 03.65. Bz, 03.65. Fd, 89.70.+C 



I. INTRODUCTION 

Quantum information must be protected against the detrimental effects of decoherence To this end 

Decoherence-Free Subspaces (DFSs) |^-[l(J have recently been proposed, alon gsid e Quantum Error Correcting Codes 
(QECCs) Q and "dynamical decoupling" and symmetrization schemes [|l5|dj8|] . A DFS is a "quiet corner" of 
the system's Hilbert space, where the evolution is decoupled from the bath and thus is entirely unitary. DFSs are a 
special class of (fully degenerate) QECCs |j| , so in order to properly distinguish between DFSs and all other QECCs 
we note that DFSs are passive codes, in that the information encoded in them may not require any active stabilization 
procedures |y^,|o). All other QECCs, in contrast, always involve an active error-detection/correction process. Exam- 
ples of DFSs have so far focused almost exclusively on the presence of a permutation symmetry of some sort in the 
system-bath coupling. The most often used example is that of "collective decoherence" where the bath 

couples in an identical fashion to all qubits, implying that all qubits undergo the same error. In this case four physical 
qubits suffice to encode a logical qubit against any collective error, and the code efficiency (number of encoded per 
physical qubits) approaches unity asymptotically Q . It was shown that the requirement of an exact symmetry can 
be lifted by allowing for a symmetry-breaking perturbation, without spoiling the DFS property significantly [p|,[To| . 
Moreover, by concatenation with an active QECC, a symmetry-broken DFS can be stabilized completely ||. While 
these results indicate that a small departure from the exact symmetry condition for the system-bath coupling is 
admissible, they leave unanswered the question of whether a DFS may exist when no assumptions are made regarding 
the spatial symmetry of this coupling. 

In this paper, the first out of two, it will be shown that under conditions which do not relate to a spatially 
symmetric system-bath coupling, DFSs may still exist. This result is exact, i.e., it is not of a perturbative nature as 
in Refs. || [Hj. Instead, it relies on the assumption that errors affecting single qubits are absent, and to lowest order 
only multiple-qubit errors are possible instead. Formally, the condition is that the qubit register is not affected by 
the full Pauli group of errors, but only by a subgroup thereof. One may then proceed to find DFSs with respect to 
this subgroup. The interesting class of system-bath interaction Hamiltonians which allow for such processes generally 
involve only multiple-qubit operators. Relevant physical systems are therefore those where the bath can only couple 
to multiple system excitations as is the case for decoherence due to dipolar coupling, e.g., in NMR p2fl . Another 
interesting class of examples are composite particles, such as bi-excitons in quantum dots/wells or Cooper pairs 
in superconductors (24). 

The structure of this paper is as follows. In section |n] we briefly review the structure of Hamiltonians pertinent to 
systems that may function as quantum computers, coupled to a decohering environment. Using these Hamiltonians, 



we recall in section EH the derivation of the operator sum representation evolution equation for the system density 



matrix. We show in particular that for a qubit system, the evolution can be expressed entirely in terms of linear 



combinations of tensor products of Pauli matrices. We then use this in section IV to derive the DFS condition 
under the assumption that decoherence is the result of a subgroup of the Pauli group. In section ^ we illustrate 
our general analysis with some exampl es, and find decoherence-free states for a number of sub grou ps. We derive 
the dimension of these DFSs in section VI. Conclusions and a summary are presented in section |VTl| . Finally, some 
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important properties of the Pauli group are summarized in Appendix |A|, and some examples of "non-generic" DFSs 
are presented in Appendix |b[ We show in a sequel paper | 25| how to perform universal fault tolerant quantum 
computation using at most 2-body Hamiltonians on the DFSs derived here. 



II. STRUCTURE OF THE HAMILTONIAN FOR A UNIVERSAL QUANTUM COMPUTER COUPLED 

TO A BATH 

This section provides a brief review of the structure of Hamiltonians relevant for a qubit system allowing for universal 
quantum computation and coupled to a decohering bath. 

The dynamics of a quantum system S coupled to a bath B (which together form a closed system) evolves unitarily 
under the combined Hamiltonian 

H = H S <E)I B + I s <8>Hb +H /; (2.1) 

where H5, Hg and H/ are the system, bath and interaction Hamiltonians, respectively; I is the identity operator. 
Let erf denote the a th Pauli matrix, a = {0, x, y, z}, acting on qubit i. The 2x2 identity matrix is denoted af. For 
K qubits the components of H can often be written as follows: 

K K 

Hs = E E £ >? + E J ^>7 + h - c -' ( 2 - 2 ) 

z—l a— x, z i^j 

where af = {af =F iaf) /2. The first sum contains the qubit energies (ef ) and tunneling elements (ef) (2(|, and the 
second sum expresses tunneling between sites i and j. Other forms are also possible, e.g., as in an anisotropic dipolar 
medium such as in solid state NMR p^| , where one would typically encounter an Ising J^ofcrJ term. A Hamiltonian 
of the form above is sufficiently general to allow for universal quantum computation by satisfying the following two 
requirements |27|-p9|: (i) Arbitrary single-qubit operations are made possible by the presence of af , which allows for 
the implementation of a continuous SU(2) rotation in the i th qubit Hilbert space, while the af term allows for the 
introduction of an arbitrary phase-shift between the |0) and |1) states. When af and af are exponentiated, they can 
be combined, using the Lie sum and product formulae |3fj] 



lim ( 



e iaA/n e iPB/n\ _ e i(aA+/3B) 



lim ( e iA/VH e iB/V* e -iA/V* e -iB/V*) = e l A ' B \ (2.3) 

to close the Lie algebra su(2), and thus to construct any evolution in the Lie group SU(2) of all possible operations 
on a single qubit p8| ]. (ii) The second ingredient needed for universal quantum computation is the controlled-not 
(CNOT) gate, which is made possible through the ability to implement each of the (nearest neighbor) JJy af^aj +h.c. 
terms. When exponentiated, such a term yields: 



with 9 oc Jijt. For 9 = 7r/4 this is (up to a phase) the "square-root-swap" operation, which when combined with 
single-qubit rotations allows for the implementation of CNOT. Alternatively, a J^afcrJ term alone is sufficient, since 
it can be used to implement a controlled-phase-shift, as is done routinely in NMR [ pl| . It is important to emphasize 
that the universal gates construction just described is but one of many different ways to achieve universal quantum 
computation. In fact, universal gates implementing logic operations directly on physical qubits (as above) are generally 
inappropriate for the purpose of fault tolerant computation |3^] . We consider a different gate construction in the sequel 
paper p5| , operating instead on "encoded" qubits, which can be used to implement universal fault tolerant quantum 
computation. For a useful survey of different universal and fault tolerant sets of gates see Ref. |3j| . 
The bath Hamiltonian can be written as 
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n B = J2 LUkBk (2.4) 
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where, e.g., for the spin-boson Hamiltonian, Bfc = bj^bfc j2(J|, and bj., b^ are respectively creation and annihilation 
operators of bath mode k. 

Finally, the system-bath interaction Hamiltonian is 

K 

H ' = E E E3>^B«, (2.5) 

i— 1 a— -\-,—,z k 

where gf k is a coupling coefRcient. In the spin-boson model one would have B,t = bfe, B,7 = h\ and B| = h\ + b^. 
Thus of <g> B J expresses a dissipative coupling (in which energy is exchanged between system and environment), and 
of B fc corresponds to a phase damping process (in which the environment randomizes the system phases, e.g., 
through elastic collisions). 

An interesting limiting case arises when the coupling constants are independent of the qubit index: gf k = g k . This 
situation, known as "collective decoherence" , arises when there is full permutational symmetry of qubit positions, and 
implies the existence of a large DFS P,pl|. Defining collective system operators S a = X)i=i a f > one can then express 
the interaction Hamiltonian in greatly simplified form as 

H coU. = ^ S"®(£ fl gBg]. 

a=+,-,z \ k / 

A case of intermediate symmetry arises when the coupling constants are equal not over the entire qubit register but 
rather only over finite clusters j = 1..C. One can then define cluster system operators S" = X)f=i a fi where Kj is 
the number of qubits in cluster j. The interaction Hamiltonian becomes 

H / lus =E E sfafetfA^V 

j = l a=+-,z \ k / 

In this case too, DFSs can be found. The point we wish to emphasize presently is that the underlying assumption 
in cluster decoherence is that of spatial symmetry in the system-bath coupling. This is to be contrasted with the 
decoherence models studied in this paper, where DFSs will be shown to arise without the need for spatial symmetry. 
Returning to the general case, H/ can be rewritten as 

K 

H ' = E E E a i ® B ? fc ' ( 2 - 6 ) 

i— 1 a— x,y,z k 

where B* fc = and Bf fc ,B^ fc are appropriate linear combinations of B^ and B^: 

B? fc = ^(9^-g+B+) (2.8) 

The qubit-coupling term in Hg can also be expressed entirely in terms of erf, where a = x, y or z. Thus all system 
components of the Hamiltonian H can be expressed in terms of tensor products of the single qubit Pauli matrices. 



III. TIME EVOLUTION OF THE DENSITY MATRIX 



The purpose of this section is to show that the evolution of the density matrix of an open system can be expanded 
in terms of tensor products of the Pauli matrices (the Pauli group), and that this follows from the structure of the 
Hamiltonians assumed above for a qubit register. This result is obvious from a formal mathematical point of view 
(since the elements of the Pauli group of order K form a complete orthogonal set for the 2 K x 2 K matrices) |Q , so 
that the reader for whom this type of argument is satisfactory may safely skip ahead to the next section. We present 
the derivation of this result here in order to motivate the appearance of the multiple- qubit errors that are the subject 
of this paper. 

We first transform to the interaction picture |3q| defined by the system and bath Hamiltonians: 
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H -> H(t) = U SB (i)HU^(t) = H 5 ®I S + I S ®H B + H/(t) (3.1) 

where 

Usb(*) = cxp [- (H S ®I B + Is^Hs) it/%] 

= exp[-jtH s /n](g)exp [-itU B /h] = Us(t)®U B (t). 

Because the system and bath operators commute, the interaction picture interaction Hamiltonian can be written as: 

K 

H/(t) = v SB (mivl B (t) = e E £°?(*)® B &(*)> ( 3 - 2 ) 

i=X a—x^y^z k 

where 

of (t) = U s (t)af U^(t) = E *if 

Ba(t)=U fl (t)B? k U^(t) (3.3) 

(see, e.g., Ref. |55| for an explicit calculation of the Xff(t) for some examples). The system-bath density matrix is 
transformed accordingly from the Schrodinger into the interaction picture (denoted by a prime): 

Ps B (t) — > p'ssit) = XJl B (t) PsB (t)U S B(t), (3.4) 

and the full dynamics is 

P'ssit) = U(ty sjJ (0)Ut(t), (3.5) 

where 



U(t) = Texp 



/ H 7 (r)dr 
Jo 



(3.6) 



and T is the Dyson time-ordering operator (defined explicitly below). From now on we work in the interaction 
picture only, so for notational simplicity the prime is dropped from the density matrices. At t = the Schrodinger 
and interaction pictures coincide. Thus, assuming that system and bath are initially decoupled so that pss(O) = 
p(0) (g> ps{0), where p and ps are, respectively, the system and bath density matrices, the system dynamics is 
described by the reduced density matrix: 

p(0) — ► p(t) = Tr B [U(t)(p(0) ® p B (0))Ut(t)]. 

Here Tr B is the partial trace over the bath. By using a spectral decomposition for the bath, p B (0) = X^P^I^K^IiQ 



this can be rewritten in the "operator sum representation" |10 36 

p(t)=J2Mt)p(0)Al(t) (3.7) 

d 

where 

A d (t) = ^(n\V(t)\v) ; d={n,v) (3.8) 
Also, by unitarity of U, one derives the normalization condition, 



1 For a bath in thermal equilibrium, \v) would be an energy eigenstate with energy E u , and p„ = exp(—(3E v )/Z, where f3 is 
the inverse temperature and Z — Tr[exp(— /3Hs)] is the canonical partition function. 
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which guarantees preservation of the trace of p: 

Tr[p(t)] = Tr[]T A d p(0) A*] = Tr[p(0) £ A* A d ] = Tr[p(0)]. (3.10) 

d d 

The {A^}, called the Kraus operators, belong to the (Banach, or Hilbert-Schmidt) space B(H) of bounded operators 
acting on the system Hilbert space, and for K qubits are represented by 2 K x 2 K matrices, just like 
Consider now a formal Taylor expansion of the propagator: 



OO / All / „t 

U (*) = E^ T / H/(r)dT 



.I + ^^|u„(t). (3.11) 



/ dt n [ dt n -i... [ di^H/fa)-- •!!/(*„)} 
Jo Jo Jo 



n 

n=l 



The Dyson time-ordered product is defined with respect to any set of operators Oj(^) as 

T{0 1 (t 1 )---O n (t n )} = T1 (t Tl )---0 T Jt Tn ) {t Tl > t T2 > ... >t Tn 



Using Eq. (2.6) we have for the terms in the above sum: 



n K n n 

II"/' E E E®-- (*i)®B&fo 

j = l i=l Q!=a;,j/,2 k j = l j=l 

where i = {ii, 12, i n }i & = {ai,c»2, ...,»„}, and k = {k\, k%, A;^,}. The important po int to notice in this 



complicated expression is that after taking the bath matrix elements (n\ ■ ■ ■ \v) [because of Eq. (3.8)], one is left with 
all possible tensor products u tj fe) over n on ^ °f K qubits. The integration and time-ordering operation will 

not change this conclusion. Thus, using, the expansion of erf (t) in Eq. fl3.3[), after a time 0(t K ) one finds the tensor 



product <2)jLi ^i* 3 ) i' e '' qubits are involved (here aj = 0, corresponding to the identity matrix, is allowed). At 
this point the entire Pauli group Pk appears (all possible A K+1 tensor products of the 3 Pauli matrices and the 
identity matrix, and the four roots of unity {±,±i} - see Appendix ^), and one has "complete decoherence" , i.e., 
multiple-qubit errors over the entire system Hilbert space. In the usual approach to QECC one does not consider such 
high orders in time since one assumes that error correction can be done quickly enough. Instead the error-analysis 
is usually confined to time evolution to 0(t) only, which leads to "independent decoherence", i.e., single-qubit errors 
affecting only one qubit at a time.^j It is possible to use multiple-error-correcting quantum codes for 0(t n ) with 
arbitrary n, but these codes are rather unwieldy (i.e., the number of encoding qubits becomes large). In the case of 
"burst errors" (a spatially contiguous cluster of errors such as / ■ • • IX ■ ■ ■ XI ■ ■ ■ I) some particularly efficient codes 
are known |^l| . 

On the other hand, a DFS that exists by virtue of a spatially symmetric system-bath coupling, is not affected by 
this proliferation of errors, which all occur in the subspace orthogonal to the DFS M. The assumption of s pati al 



symmetry manifests itself in restrictions on the coefficients gf k appearing in the interaction Hamiltonian [Eq. (2.5)]. 
For example, as mentioned above, collective decoherence corresponds to the condition gf k = g k Vi, i.e., the bath 
cannot distinguish between the qubits ||. In this paper no such spatial symmetry assumptions will be made. Instead, 
only multiple-qubit errors will be allowed to lowest order instead of single-qubit errors. This condition will be defined 
more precisely in the next section. 

As for the Kraus operators, it can be seen from the calculations above that they may be expanded as sums over 
tensor products of the Pauli matrices: 

Ad(t) = a dAt)Pn (3.12) 



2 See, however, Ref. J38| for a discussion of Kraus operators represented by non-square matrices. 

3 In fact spatially correlated errors can also be dealt with by QECCs wUt. 
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where p„ € Pk- The Kraus operators thus belong to the group algebra (the space of linear combinations of group 
elements) of Pk if^j - As alluded to in the beginning of this section, that this expansion is possible actually follows 
simply from the fact that the Pauli group forms a complete orthogonal set (with respect to the trace inner product) 
for the expansion (with complex coefficients) of arbitrary 2 K x 2 K matrices. However, here we have seen how the 
expansion in terms of the Pauli group (rather than some other basis) is physically motivated by virtue of the structure 
of the Hamiltonian. 

A simple example will now serve to illustrate the point made above about multiple-qubit errors. Consider an 
interaction Hamiltonian of the form Hj — J2i=i a t ® ( on two qubits). Some algebra suffices to show that then 
Ad{t) — c (t)Is + ci(t)erf + C2(i)erf + c\iif)a\ ® erf. In this case the single-qubit errors erf, erf appear, as well as 
the multiple-qubit error o~\ <§5 erf . This situation does not allow for the appearance of DFSs (unless spatial symmetry 
is present). Alternatively, consider the interaction Hamiltonian H/ = {o\ <g> erf) ® _B 12 + (erf ® erf) ® -B34 (on four 
qubits). In this case one finds Ad{t) = c (i)Is + Ci2{t)af (£> erf + C34(t)crf £g> erf + Ci234(t)af <g> erf <g) erf ig) cr|. Thus 
only multiple-qubit terms appear, and as will be shown below, this allows for the existence of non-trivial DFSs, even 
though no spatial symmetry assumptions were made. 

An important example of this correlated type of system-bath interaction is the dipolar-coupling Hamiltonian, 
relevant, e.g., to decoherence resulting from spin-rotation coupling in NMR |^2|J^| The dipolar Hamiltonian for a 
system of spins interacting with a bath of rotations is: 

Hi = E 7? ' - 3 ' *>) ' r ^)] > ( 3 - 13 ) 

j,k i k 



where jj is the gyromagnetic ration of spin j, rjk is the distance beween spins j and fc, and er is the vector of Pauli 

a/3 

9jk 



matrices. Introducing an anistropy tensor g°*f~ , this can be rewritten as: 



1 



E at^^o-^Y-^, (3.14) 



j,k k a,/3=-l 

where are the spherical harmonics, and Eventhough only multiple-qubit terms appear here it is necessary 

to further impose anisotropy in order to obtain an example with a non-trivial DFS, as we discuss in more detail in 



Sec. V A4 . This is the case, e.g., when only crj <g> a\ terms remain (i.e., = S a oo~/3ogjk), coupled to Y® rotations 



With these observations, we are now ready to study the question of DFSs in open systems without spatial symmetry 
in the system-bath couplings. 

IV. DECOHERENCE-FREE SUBSPACES FROM SUBGROUPS OF THE PAULI GROUP 

We begin this section by recalling the condition for DFSs within the framework of the Kraus operator-sum repre- 
sentation, derived in Ref. || . We then analyze the conditions for the appearance of DFSs when the errors are spanned 
by a subgroup of the Pauli group. The result is summarized by a theorem presented at the end of the section. 

A. Condition for Decoherence-Free Subspaces 

A DFS is a subspace H — Span{|j)} of the full system Hilbert space Hk over which the evolution of the density 
matrix is unitary. Necessary and sufficient conditions for a DFS were derived in the Markovian case in Ref. Q and 
in the exact (non-Markovian) case in Ref. [|| . A formulation of the exact DFS condition was given in terms of the 
operator sum representation in Ref. and will be briefly reviewed. 

Let { I j) } be a set of system states satisfying: 

A d \j) - c d fj\j) Vd, (4.1) 

where U is an arbitrary, ^-independent but possibly time-dependent unitary transformation, and a a complex con- 
stant. Under this condition, an initially pure state belonging to Span[{|j)}], 



We thank Prof. Dieter Suter for suggesting this example. 



G 



\4>in) = 

j 



will be decoherence-free, since: 



\4> d ) = A d |^in) = jjCdUlj) = CdU|^i E 



so 



Pout 



d 



CjUl^mXV^lt^cS = U|Vin)(V'in|U t 



where we used the normalization of the Kraus operators [Eq. ( pTg| )] to set 5Z)rf I Cd 1 2 = 1- This means that the 
time-evolved state p out is pure, and its evolution is governed by U. This argument is easily generalized to an initial 
mixed state p- m = Y]jf Pjj'\j)(j'U m which case p ut = Upi n U^. The unitary transformation U is a "gauge freedom" 
which can be exploited in choosing a driving system Hamiltonian which implements a useful evolution on the DFS. 
In the interaction picture used in the previous section, U can be made to disappear by redefining all Kraus operators 
as U^Arf. The calculation above shows that Eq. (4.1) is a sufficient condition for a DFS. It follows from the results 
of RefSjJ^,^ that it is also a necessary condition for a DFS (under "generic" conditions - to be explained below). 

Eq. (4.1), however, does not seem to be a very useful characterization of a DFS if one does not know the explicit form 
of the Kraus operators (in general this cannot be found in closed analytical form, althoug h th ey can be determined 
experimentally jjjj). When the Kraus operators derive from a Hamiltonian, as in Eq. (3.8), an equivalent DFS 
condition is || : 



S a \j) = a a \j) 



Va, 



(4.2) 



where the system-bath interaction Hamiltonian is written as Hj = ^2 a S a (g> B a [compare to Eq. (2.6)], with {S Q } 
being the system operators. To make use of this last DFS condition, one needs to introduce assumptions about the 
structure of system- bath coupling, and this is how one is led to spatial symmetry considerations Here, however, 
the DFS condition of Eq. ([O]) will be considered directly, based purely on the expansion of the Kraus operators in 
terms of the Pauli group elements, and without resorting to an explicit form for these operators. 



B. Representation Theory Construction of Decoherence-Free States 

When the Kraus operators are viewed as operators in the algebra of the Pauli group, the DFS condition [Eq. (|4.l|)] 
has a natural interpretation: the decoherence free states {\j}} belong to the one-dimensional irreps of the Pauli group. 
Motivated by this observation we now consider a group representation theory construction of decoherence-free states. 
The general criterion for the reducibility of a representation {r(G n )}^ =1 of a finite group Q = {G n } of order N is 




N 

£| X [r(G„)]| 2 >iV, (4.3) 

n=l 

where x IS the character of the representation T [trace of the matrix T(G n )]. If equality holds, then the representation 
is irreducible. 

The full Pauli group Pg is irreducible over the Hilbert space Hk of K qubits: since all Pauli matrices are traceless, 
only the four elements proportional to the identity matrix contribute (see also Appendix |a|) : 

4 A'+1 

\xlPn}\ 2 = \2 K \ 2 + | - 2*| 2 + K2 A '| 2 + | - ,2*| 2 = 4*+\ 

n=l 

which is just the order of Pk (generally the direct product representation of irreps of any direct product group is 
itself an irrep of that group [Q). 

Now we come to the central assumption setting the stage for the DFSs considered in this paper: what if the Kraus 
operators belong to the group algebra of a subgroup Q of Pk? The motivation for this situation could be: (i) The 
case in which either only higher order errors occur, such that first-order terms of the form I®---®I®a"®It&---tg)I 
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are absent in the Pauli-group expansion of the Kraus operators, or (ii) only errors of one kind, either a x , <r v , or a z 
take place. Case (i) would imply that either: (a) There are certain cancell ations involving bath matrix element terms 



such that first-order system operators are absent in the expansion of Eq. (3.11). This would be a rather non-generic 
situation, involving a very special "friendly" bath; or (b) The system-bath Hamiltonian is in fact not of the form 
in Eq. ( |2.6| ), but rather involves only second order terms such as af ® (identity on all the rest). ^| Case (ii) is 
applicable in, e.g., the case of pure phase damping (relevant to NMR and optical lattices using cold controlled 
collisions p5[|), where o~ z errors are dominant. 

In the subgroup case under consideration, we may find non-trivial irreducible representations (irreps) of Q over Tlx 
(a so-called "subduced" representation [^U). This situation can be interesting especially if there exist 1-dimensional 
irreps, as known from the general theory of DFSs As will be shown next, the recipe for finding these DFSs uses 

the standard projection operators from elementary group representation theory. The projection is onto the subspace 
transforming according to a particular irrep. 

First, recall the multiplicity formula for unitary irreps (which we can always assume in this case since the Pauli 
group is finite): 



N 

1 

TOfc 



1 N 

-]Tx[r fc (G„)]*x[r(G„)], (4.4) 



where rrik is the number of times irrep T k appears in the given reducible representation; \ [r fc (G n )] is the character 
of the r fc irrep on the group element G„; and x [r(G ra )] is the character of G„ in the given reducible representation T. 

We denote a set of (orthonormal) basis-states transforming according to an irrep r fc by {IV'i), ■ ■ ■ > l^d fc )}- These 
states span the invariant subspace of the irrep T k and transform according to: 

G n \tf) =J2T k (G n ) vfl \ip k u ). (4.5) 
i/=i 

Furthermore they obey the orthogonality relation: 

(i>l\i> k v ) = 5 lk 6^ (4.6) 

Next, a projection operator onto the subspace belonging to the d k -dimensional irrep k is given by the appropriate 
sum over group elements f4~i|| : 

d N 

P l = 4Y1 r k (G n )*G n ; n,u=l, d k (4.7) 



N ^ 

n=l 



and has the following properties: 



PlM)=hi5vx\^t). (4.8) 

To obtain a set of (orthonormal) basis states {|V , i)i • • • > IV^ fc )} transforming as a set of partners in the basis for T k from 
an arbitrary state \<j>), one can apply the set of operators {P k v } for a fixed v (such that P k v \4>) ^ 0) and renormalizc 
the states thus obtained. Every state \<f) can be expanded in terms of basis states for the constituting irreps T k as 

(4.9) 

k u=l 

where P k v \(j>) = k \il) k ) and the summation over k is over inequivalent irreps p4j . 

Let us now consider the effect of applying the operators = J2 n o-d,nG n from the group algebra to an arbitrary 
state \6): 



5 Note that in this case the expansion of the Kraus operators in terms of tensor products of Pauli matrices, Eq. (3.12) remains 
valid. 
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dk N d k d k 

E E e iM^D = E E E " E r^wi^) 

k jU=l n— 1 fe /x=l z^— 1 



(4.10) 



We would like to find the conditions such that this transforms into the DFS condition, Eq. (4.1). Consider the case 
when r fc are all 1- dimensional irreps, possibly appearing with multiplicity m^: 



In this case the indices /x, v are irrelevant and we will omit them. Then: 

N 



A d |0) = E arf .™E^ fc i^ fe ) 



(4.11) 



(4.12) 



For \4>) to be a decoherence free state, one would like to have this proportional to \<fr) = J^k S k \ij] k ) [as in the original 
expansion of Eq. ([0])]. However, this does not work because of the presence of j k in the sum. We thus see that the 
initial function \(j>) must be restricted to be one of the basis-states \tp ). Then, with 



N 

E 



a d,nJ„ 



(4.13) 



we have finally: 



A d \^) = c k \^ k 



(4.14) 



At this point it is useful to introduce another index z for the multiplicity of the irrep k, i.e. z = 1 . . . mfe. The Hilbert 
space of -fT-qubit states splits into invariant one-dimensional subspaces V k that are spanned by (fixed) basis states 
\ip k ). Each of the \ip k ) in Eq. (|4.9D is a linear combination of the \i/j k ): 



iv>*>=E^>' 



(4.15) 



[Because of Eq. (^9|), the Q k z depend on the initial state \(f>}.] Thus for \(f>) to be a decoherence-free state, it is allowed 
to be an arbitrary superposition inside copies of a given irrep (different z's), but not to be a superposition between 
different irreps (different fc's). In particular we have within each copy of the irrep T k : 



A d \ip k ) = c k M 



k ) ■ 



Z = L.TOfe. 



(4.16) 



This is just the DFS condition, Eq. ( p~i"| ) with the {1^)} being the basis states for the DFS. Therefore Eq. ( |4.1l| ) 
is a sufficient condition for a DFS, provided that our initial state satisfies the condition that it is a superposition of 
states within a fixed irrep, Eq. ( 4.15|). 

It will now be shown that Eq. ( [4.11 ) is also a necessary condition for a DFS under the "genericity" assumption 
that the error coefficients {ad,n} are arbitrary. In other words, it will be shown under these conditions that, if a set 



of basis-states {\j}} satisfies the DFS-condition Eq. ( |4.l| ) then the belong to the invariant subspace of some 

one-dimensional irrep of our subgroup. 



Assume that the have been redefined to incorporat e th e (constant) unitary transformation U such that Eq. (4.1) 
becomes A^j) = Cd\j)- Expand the state \j) as in Eq. (4.9): [] 



b> = EE^, 



(4.17) 



k M=l 



where P^\4>) = 0£ fe |V>£). Now using Eq. ( fiic|) 



For notational simplicity we avoid introd ucing another index for the multiplicity of the irrep here. That such superpositions 
are allowed for DF states is clear from Eq. (4.15). 







dk 

A d \~ 3 ) = c d \]) = c d j2J2 d l k \0 ( 4 - 18 ) 

k fi=l 

= e E ^ A <*i^> = E E E a ^ G «i^> 

k k n 

= E E 9 l k E E r fc (G„)A M i^) (4.19) 

n A=l 



and taking inner products [using (4.6)] 



(^iA d |j) = qE E d l k & l M) = 

k 

= EE^E^E rl ( G «w^i^) = E^E a ^ r '( G «w- (4-20) 



Using this result we would like to show that the T l (G n ) that appear here must be one-dimensional irreps. Let us 
establi sh "g eneric" conditions for this purpose. 



Eq. ( 4.20 ) can be rewritten as an eigenvalue equation: 

4$ = c 4< ( 4 - 21 ) 

where 

A l d = J2 a cLnT l (G n ) (4.22) 

n 

? = (4.23) 



The vector #J may be zero for a given irrep r', in which cas e Eq. ( 4.21 ) is trivially satisfied. Let us assume t his i s 



not the case for some I [it cannot be the case for all I, by Eq. (4.17)]. Then the most general way in which Eq. (4.21) 



-3 



can be satisfied, is for - to be an eigenvector of A d for all codewords with eigenvalue c d . However, while this is 
the most general condition, it is non-generic. By generic we mean that we take the errors to be arbitrary, i.e., we do 
not want to make any assumptions on the a d ^ n . Now, if the eigenvalue eqnarray were to be satisfied, the vector of 

coefficients would have to be "special". In other words, it would have to be adjusted to be an eigenvector of A d . 
To make this adjustment would require two conditions: (i) Having a priori knowledge of the a d ^ n , (ii) Being able to 

control ~9j . We would like to avoid assuming (i) because fine-tuning the bath is physically unacceptable. In contrast, 

control of is certainly desirable. However, we would like to avoid the situation where only certain special choices 

o t ?, — e * sp ec ifl c b a tk P _, yl e!d decode free flfl We eo„c lu de to „o,d 

fine-tuning of the bath parameters and/or special initial conditions, A l d must be proportional to the identity. But 
since T l is an irrep this is only possible if it is one-dimensional, i.e. T l (G n )^ v = 7^, fx, v = 1 and c d — J2 n a d,nl l n - In 
addition we see that c d can only be Z-independent if the DFS states \j) are spanned only by basis states of copies of 
the same irrep T . Q.E.D. 
To summarize: 

Theorem 1: Suppose that the Kraus operators belong to the group algebra of some group Q = {G n }, i.e., = 
Tlin=i a d,nGn- If a set of states {\j}} belong to a given one- dimensional irrep T k of Q, then the DFS condition 
A d \j) — c d \j) holds. If no assumptions are made on the bath coefficients {a d . n }, then the DFS condition A d \j) = c d \j) 
implies that \j) belongs to a one- dimensional irrep T k of Q. 

For completeness we give in Appendix [b| an example of the "non-generic DFSs" , which result from "accidentally" 



satisfying Eq. (4.21) with irreps of dimension greater than one. 



7 This statement of what are generic conditions that lead to a DFS is very similar to that in Ref. 
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V. EXAMPLES OF SUBGROUPS WITH DECOHERENCE FREE STATES 



The general considerations from the previous section will now be illustrated with some examples. To simplify the 
notation, let X,Y,Z represent the a x ,<r y ,a z Pauli matrices, and let us drop the tensor product symbol (i.e., let 
ZI = Z (g) I, X 2 = X (g) X, etc.). Also, we will ignore normalization factors in this section. 



A. Abelian Subgroups 

The simplest non-trivial example of a subgroup is found already for K = 2 qubits: 

Qz = {I 2 ,ZI,IZ,Z 2 }. (5.1) 

This subgroup (generated by ZI and IZ) describes phase damping. 

As another simple example, let K = 4 qubits and consider the following subgroup: 

Qx = {I 4 ,X 2 I 2 ,I 2 X 2 ,X 4 }. 

Physically, this would correspond to the error process where bit flips happen on certain clusters of two or four qubits 
only (note that XIX I and IX IX were left out - this case will be considered in the sequel paper |j25f). 
Another example is 

Q 4 = {I 4 ,X 4 ,Y 4 ,Z 4 }, 

with all Pauli errors occurring just on clusters of 4 qubits. Qz, Qx and Q4 are isomorphic and Abelian. All elements 
of these subgroups, except I 4 , are traceless. I 4 has trace 16, so that J2n=i lx[r(G„)]| 2 = 256 > 4 and thus the natural 
representation of these subgroups on 4 qubits is reducible. Since they are Abelian, they have only 1-dimensional 
irreps. These irreps are given in the following table, expressed in terms of the elements of Qx- 





/ 


XI 




X 


r 1 


1 


1 


1 


1 


T' 2 


1 


1 


-1 


-1 




1 


-1 


1 


-1 


r 4 


1 


-1 


-1 


1 



Motivated by Theorem 1, this reducibility implies the existence of DFSs, as long as the Kraus operators belong to 
the group algebra of these subgroups. 



1. The Subgroup 



Consider the case of Qx, i-c, assume that the Kraus operators can be written as 

A d = a dfi I 4 + a dA X 2 I 2 + a d . 2 I 2 X 2 + a d , 3 X 4 (5.3) 

[the coefficients a d j are of course const raine d by the normalization con dition Eq. ( |3.9| )]. 

Using the general arguments of Sec. [VB and in particular Eq. ( |4.7| ), we can just read off the matrix elements of 
the 4 (1-dimensional) irreps from the table in Eq. (5.2). Thus the 4 projection operators are: 



P 1 = I 4 + X 2 ! 2 



P 3 = I 4 



X 2 ! 2 



vx 2 

j2 x 2 



X 4 
X 4 



P 2 



I 4 + X 2 ! 2 - I 2 X 2 - X 4 



P 4 = I 4 



X 2 ! 2 



I 2 X 2 



X 



(5.4) 



The multiplicity of each of the four 1-dimcnsional irreps in t he r educible representation generated here by the K = 4 
qubits, is 4. To see this, recall the multiplicity formula Eq. (4^4). In the present case, the given representation yields 



X = {16,0,0,0} (for I 4 ,X 2 I 2 ,I 2 X 2 ,X 4 respectively) and so, with X k (I 4 ) = l,m k = jX k (I 4 )^ = 4 for all k. 

Now, let us explicitly find the decoherence free states. To do so we can pick an arbitrary, convenient 4-qubit state 
and project it onto a given irrep. For example, starting with 1 0000) : 
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P 1 10000) = |0000) + |1100) 
P 2 |0000) = |0000) + |1100) 
P 3 |0000) = |0000) - |1100) 



loon) + |iiii) = ivi) 

|0011) - |1111) = i^f> 
|ooii) -|ini) = |v?) 

P 4 |0000) = |0000) - 1 1100) - |0011) + |1111) = IV'i). 



(5.5) 



Each of these 4 states belongs to a different irrep, and thus to a different DFS, which can be verified by applying an 



arbitrary Kraus operator, as in Eq. (5.2). E.g., 



A d |^i) = <z d , Q (10000) + 
+o dl i(|1100) 
+a d , 2 (|0011) 
+a d , 3 (|llll) 
= {a>d,o + a dj i + 



|1100) + |0011) + |1111>) 
f |0000) + 1 1111) + |0011)) 
f 1 1111) -h |0000) + |1100» 
f |0011) + |1100) + |0000)) 
ad,2 + a d ,3) Vl>\)- 



(5.6) 



Similarly: 



(a d ,o + a d ,i - a d , 2 - a d ^) \tjjj), 
(ad,o - a d ,i + a d ,2 - o-dfl) |Vi)) 
(a d ,o - 0>d,i - a d,2 + a d ^) \tpf). 



(5.7) 



This is in agreement with Eq. (4.16) 



Now, recall that each irrep appears 4 times. This means we should be able to find 3 more independent states 
belonging to each of the irreps. Indeed, by performing projections on the states |0001), |0100), 1 1001) (using |0010) 
and 1 1000) does not produce new states) we obtain the complete basis for the DFSs. E.g., 



and again 



AdlVz) 



P 1 10001) = |0001) + |1101) + |0010) + |1110) = \ipl), 
P^OIOO) = |0100) + 1 1000) + |0111) + |1011) = \ipl), 

p^iooi) = I looi) + |oioi) + |ioio) + |oiio) = IV'I), 



[a dfl I 4 + a ±1 X 2 I 2 + a d . 2 I 2 X 2 + a d . 3 X 4 ] (|0001) + |1101) 
(a-dfi + a d ,i + a da + a d , 3 ) \ipl), 



|0010) + 1 1110)) 



(5.8) 



(5.9) 



with similar results for the other states. All of this is in agreement with the general results of Sec. [V B| . Finally, we 
may consider an arbitrary superposition of decoherence free states taken from the multiple appearances of a given 
irrep, \cj) k ) = ^zl^z)' an< ^ tn ' s wm a g am be decoherence- free. 



2. The Subgroup 
In this case the Kraus operators can be written as 



A d = a d . I 4 + a dA X 4 + a da Y 4 + a d , 3 Z 4 . (5.10) 



Again, using the general arguments of Sec. [IV B| , in the case of Q4 we can just read off the matrix elements of the 
4 (1-dimcnsional) irreps from the table in Eq. ( |5.2|) . Thus the 4 projection operators are: 

P 1 = I 4 + X 4 + Y 4 + Z 4 P 2 =I 4 + X 4 -Y 4 - Z 4 

P 3 = I 4 -X 4 + Y 4 -Z 4 P 4 = I 4 -X 4 -Y 4 + Z 4 . (5.11) 



Using the multiplicity formula, Eq. (4.4), the given representation again yields \ — {16,0,0,0} (for J , X 4 ,Y 4 , Z A 
respectively) and so once more mjt — |x fe (/ 4 )16 = 4 for all k. 

To find the decoherence free states let us start again with 1 0000) . We find: 
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P 1 1 0000) = 2(|0000) + |1111» = |Vi) 

P 2 |0000) = |0000) + |1111) - |1111) - |0000) = 

P 3 |0000) = |0000) - |1111) + |1111) - |0000) = 

P 4 |0000) = 2(|0000) - |1111» = |Vi)- (5.12) 

The vanishing of the projections of P 2 and P 3 implies that |0000) has no components in the irrcps T 2 and T 3 . Thus 
a different starting state is needed, e.g., |0001). Then 

P 2 |0001) = 2(|0001) + 1 1110)) ee \iPl) 

P 3 |0001) = 2(|0001) - 1 1110)) ee |i/> 3 ). (5.13) 

That these states are decoherence free, is again easily verified by application of an arbitrary Kraus operator, e.g.: 

AalVi 2 ) - [a d , I* + a dA X 4 + a d . 2 Y 4 + a d , 3 Z 4 ] 2 (|0001) + |1H0)) 

= (o<i,o + a d.i - a d ,2 - a d $) IV'i), (5-14) 

etc. The full DFS corresponding to the projection P 1 is found by applying P 1 to the initial states |0011), |0101), 1 1001) : 

P^OOll) = 2(|0011) + |1100)) ee |v4) 
P^OlOl) = 2(|0101) + |1010» ee |V>a) 

P^lOOl) = 2(|1001) + |0110)) ee |Vl), (5.15) 

in addition to above. 

Since the decoherence process described by Q4 is different from that of Qx, the decoherence free states are, not 
surprisingly, different in the two cases. 



3. The Subgroup Qz 

As another example of an Abelian subgroup, assume now that the Kraus operators, for K = 2 qubits, can be 
written as 

A d = a dfi I 2 + a dA ZI + a d . 2 IZ + a d ^Z 2 . (5.16) 
The 4 projection operators are thus: 

P 1 = I 2 + ZI + IZ + Z 2 P 2 = I 2 + ZI - IZ - Z 2 

P 3 = I 2 - ZI + IZ - Z 2 P 4 = I 2 - ZI - IZ + Z 2 (5.17) 

In this case, the given representation on 2 qubits yields \ = {4,0,0,0} (for I 2 ,ZI,IZ,Z 2 respectively) and so 
TOfc = i% fc (7 2 )4 = 1 for all k. Thus as expected (since the representation is 4-dimensional) , the multiplicity of each of 
the four 1-dimcnsional irreps is 1. 

Let us again explicitly find the decoherence free states: 



P^OO) 


= 4|00) EE 


l^ 1 ) 


P 2 |01) 


= 4|01) = 


IV 2 ) 


P 3 |10) 


= 4|10) = 


IV 3 ) 


P 4 I11) 


= 4|11) = 


IV 4 ) 



And indeed: 

A d |^ fe ) = (a d:0 + ad.i + a da + a d>3 ) \ip k ), k = 1, 4. 

This means that each of the 4 "computational basis states" \ip k ) is by itself a DFS. However, since these DFSs belong 
to different irrcps, a superposition is not decohcrence-free. This agrees with the well known fact that phase damping 
leads to decay of the off diagonal elements of the density matrix in the computational basis, but does not cause any 
population decay. 
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4- The Subgroup Q^z 



As a final example of an Abelian subgroup, let us return to the anistropic dipolar-coupling Hamiltonian [Eq. ( |3.14 )] 



discussed in Sec. [II. Note first that it is necessary to transform from the a basis used there to a x ' v in order for our 
Pauli group-based discussion to apply. Having done that, it is clear that unless anisotropy is imposed this Hamiltonian 
generates the entire Pauli group, since all bilinear combinations a a ® appear in it. Assume therefore that we have 



a 4-spin molecule constrained to rotate only about the z-axis. This amounts to setting g^P = 5 a o5po9jk i n Eq- (3-14) 



so that only crj <g) ct^ terms remain. The corresponding subgroup is 

Q-iz = {IHI, ZZII, ZIIZ, IIZZ, ZIZI, IZZI, IZIZ, ZZZZ}. (5.19) 

To find the DFS under Q 2 z, construct the projector P 1 — ^J2 q eQ 2Z 1 corres P on ding to the identity irrep of Q 2 z- 
Applying this projector to the initial states |0000) and 1 1 1 1 1) we find a 2-dimensional DFS, spanned by these two 
states. This DFS thus encodes a single qubit. 

B. Non-Abelian Subgroups? 

It would have been interesting to find examples of non- Abelian subgroups which have 1-dimensional irreps and thus 
support a DFS. However, no such subgroups exist in the case of the Pauli group, as we now prove. 

Each two elements of the Pauli group Pr- either commute or anticommute (Appendix |a|) . Let Q be a non- Abelian 
subgroup of Pk- Then there must be at least two elements of Q, say q\ and q 2 , that anticommute. Assume that 
the state \i) belongs to a 1-dimensional irrep T of Q. Then T(qi)\i) — ci\i) and T(q 2 )\i) = C2 1 i) , where c%,C2 
are numbers. Now, by assumption Tfaqi) = T(—q±q2). Therefore T(qiq 2 )\i) = T (qi)T (q 2 )\i) — cxC2\i), and also 
f(?ig2)|i) = f(-q 2 qi)\i) = f(-q 2 )T(q 1 )\i) = cif (-g 2 )|i). If f(-g 2 ) = -f(g 2 ) then we have f(qiq 2 )\i) = -cic 2 \i) 
so that cic 2 = — cic 2 . This implies that at least one of c\ and c 2 is zero. However, this cannot be true since the 
representation is unitary. Is there another possibility? Note that r(— q%) = T(— Iq 2 ) = T(— I)r(q 2 ), so the question 
boils down to the value of a in f (-1) = af (I). But since (-I)(-I) = I it follows that f (-l)f (-1) = f (I) = 1, so 



that r(— I) = ±1. Assume then that the other case, T(— I) = 1, holds. Let us use Eq. (4.4 ) while recalling that only 
the four multiples of the identity have non-vanishing trace: 



m k = 



1 N 

-^x[r fc (p„)]*x[r( P „)] 



N 

n=l 



= I ( x [r k (i)Y (2 K ) + x [r k (-i)Y (-2 K ) + x [r k (ii)]* (i2 K ) + x [r*Hi)]* (-i2 K )) . (5.20) 

Since the irrep T k is 1-dimensional, x \^ k ] = T* 1 , i.e., the character is the element itself. Now let r fe = T. Then since 
r(— I) = 1, and using T k (—iT) — T(— I)r(il), we find fh = 0. Therefore such irreps do not appear at all. 

Thus an anti-commuting pair of elements in Q is incompatible with a 1-dimensional irrep, so that if Q has a 
1-dimensional irrep, it must be Abelianj^] 

Recall that the DFS condition of theorem 1 applies to arbitrary groups. Groups other than the Pauli group may 
support non- Abelian subgroups with 1-dimensional irreps (the above proof relied strongly on a property specific to the 
Pauli group, that its elem ents cither commute or anticommute). However, at least within the Hamiltonian framework 
expounded in Sees. || and III, it is the Pauli group which appears naturally for the group algebra to which the Kraus 



operators belong. 

VI. DIMENSION OF THE DECOHERENCE-FREE SUBSPACES 

We showed in the previous section that for the Pauli group, DFSs can exist only for Abelian subgroups. This 



observation allows us to calculate the dimension of these DFSs. Recall from the general discussion in Sec. IV B that 



in the generic case a superposition of states belonging to different irreps will decohere, whereas a superposition of 



8 We thank Dr. P. Zanardi for discussions regarding this point. 
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states within copies of a given irrep will be decoherence-free (see also the examples in the previous section). Also, by 
the Abelian property, each such copy only supports a single decoherence free state. Hence the dimension of the DFS 
associated with a given irrep T k is simply its multiplicity rrik- 

Let Q be an order- N Abelian subgroup of the Pauli group on K qubits. Using Eq. ( |5.20| ) and r fe (-I) = ±1 again, 
we have two (and only two) cases: (i) If T k {— I) = 1 then = 0, so such irreps do not support a DFS. (ii) If 
T k {-I) = -1 then 

m k = 2 K+2 /N. 

This shows that all irreps that support a DFS have the same multiplicity, and thus all these DFSs have the same 
dimension. 

If the subgroup doe s not include elements with the ±1, ±i factors, as in the examples in Sec. [v|, then only the term 
r fe (I) appears in Eq. (5.20), and consequently 

m k = 2 K /N no {±l,±i} factors. 

In any case, the dimension of the DFS is inversely proportional to the order of the subgroup. This implies a trade-off 
between the number of errors that can be dealt with by the code (N) and the number of decoherence free qubits 
(log 2 m k ). 

As an interesting corollary we see that largest Abelian subgroup of the Pauli group has order 2 K+2 (since > 1 
implies TV < 2 K+2 ). Examples of such subgroups are: 

• The group generated by all the single qubit X's (or Y's or Z 's), with ±l,±i. 

• The group generated by XXII. ..II, YYII...II, ZZII...II, IIXXII..JI, IIYYII...II, IIZZII...II, 
II...IIZZ, with ±l,±i 

These groups support only 1-dimensional DFSs. The last group is relevant for errors due to exchange on pairs of 
identical qubits |4q1, and we see that the corresponding decoherence free state is automatically immune to exchange 
errors. (See Ref. ]2l| for a discussion of protection of DFSs against exchange errors arising in the spatially correlated 
collective decoherence case.) 



VII. SUMMARY AND CONCLUSIONS 



Decoherence-free subspaces (DFSs) are associated most commonly with the existence of a spatial symmetry in 
the system-bath coupling, as in the collective decoherence model. Here we have considered the case when no such 
symmetry is assumed, and have shown that one can nevertheless find DFSs under certain conditions. The essential 
assumptions are that either to lowest order only multiple-qabit errors are possible, meaning that the bath can only 
couple to multiple system excitations; or, that only one type of error process (such as phase-damping) occurs, which can 
be relevant for the NMR quantum computer schemes and optical lattices (or any other realization where scattering- 
induced phase-shifts are the dominant decoherence mechanism). In either case, instead of the full Pauli group of 
errors, only a subgroup needs to be considered. Barring certain non-generic cases, the DFSs then correspond to states 
that transform according to the 1 -dimensional irreducible representations of such a subgroup. This characterization 
of DFSs, while formally similar to previous results, is different in that it trades the assumption of spatial symmetry 
for one of multiple-qubit coupling to the bath. 

We show in a sequel paper p5[ ] how to perform universal fault tolerant quantum computation on the DFSs found 
in this paper using only one- and two-body Hamiltonians. It would further be desirable to identify in detail the 
physical conditions under which the Pauli subgroup model is relevant for current proposals for quantum computers. 
An important example we have discussed is the dipolar-coupling induced decoherence in NMR. 
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APPENDIX A: THE PAULI GROUP 



The Pauli matrices are: 



a ^ I= {l°i) a * = (ioj a y = (°i oj CT * = (o-°i)- ^ 

They have the following properties: 

<?l=I a = 0,x,y,z 

O a O{i — ieap-fG-f + S a fjl. 

Tr(<r a )=0 a = x,y,z (A2) 

The Pauli group of order K is the set of all 4 K+1 possible tensor products of K of the Pauli matrices and ±, ±i: 

( k \ 



P K = ±, ±i < Q$) <r a , k > • (A3) 

U=l J a 

Some of its useful properties are 

• Let pi,P2 € Pk- Since either [<7 a ,fc, <7/3,k] — or {<7 Q ,fc, 073^} = it follows that 

either \pi,p 2 ] = or {pi,p 2 } = 0. (A4) 

• Since a a are all unitary, so are all p G Pk- 

• Since a a are all Hermitian but we allow for ±i factors, p e Pk is either Hermitian or anti-Hermitian. Thus if 
pE P K then pt e p^. 

• Since Tr(A ® 5) = TrA x Tri3, the only elements in Pr- which are not traceless are the four ±, ±i multiples of 
the identity, and each has trace 2 K . 

APPENDIX B: EXAMPLES OF NON-GENERIC DECOHERENCE-FREE SUBSPACES 

We will show here an example of a DFS that arises out of a two-dimensional irrep of a non-Abelian subgroup, in 
the "non-generic" case. 

Let us consider the non-Abelian 8-element subgroup Q 8 = {±III,±XXI,±IZZ,±iXYZ}. In this standard 
representation it is reducible and splits into 4 copies of a two-dimensional irreducible representation of Q$. Since 
there is just one irrep, we drop the irrep-index k on T k etc. The two-dimensional representation of Qs is the following: 

r(±///) = ±(j J) r(±xx/) = ±(; I 

T{±IZZ) =±(J_° 1 ) T{±iXYZ) = ± ( J - 1 ) . (Bl) 

The 8-dimensional Hilbert space is split into 4 irreducible subspaces V z (corresponding to the 4 copies of T) spanned 
by 

^ 1 = (|^),IV'i 1 )) = (|ooo),|iio)) 
^ 2 = (IV-o 2 )J^)) = (|ni)Jooi)) 

V 3 = (|V>o^i>) = (|100),|010)) 

F 4 = (|Vo 4 ),l^)) = (|011),|101)). (B2) 
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On each of these two-dimensional subspaces the group acts like T. A codeword in the DFS can be expanded as 
\j) = J2t=i X^ t =o ^z^JV'p-i Let us take as our code just the first basis-vector of each irreducible subspace, i.e., 

C = {|1>, |2), |3>, |4)} = W a ) : z = 1 . ..4} = {|000>, |111), |100), |011>}. (B3) 



Denoting the vector of coefficients as 9 3 z = (Q z , Gr z -jj, this means that 91 = (1,0) and 9 J f z = (0,0). In this case we 
can show that there are Kraus operators that satisfy the DFS-condition on the code, by searching for matrices 

Ad which have 9° z as eigenvectors. An example is: 

M=(^t\ A 2 =( C * d :), (B4) 



ei J A \ e 2 _ 

with the conditions c\d\ + c^di = 0, |ci| 2 + | C2 1 2 = 1 and |di| 2 + |g?2| 2 + | ex | 2 + | e 2 1 2 = 1 for normalization [Eq. 
( refeq:OSRnorm)]. The corresponding Kraus-operators are: 

= £L±£i m + ?± XXI + 2l^l IZZ + <± lXYZ (B5) 
2 2 2 2 v ' 

Co + do Co — eo do 

A 2 = 2 2 III + —XXI + — -IZZ + —iXYZ. 

2 2 2 2 2 

The code C is a DFS. It is the particular equality (i.e., the "conspiring", non-generic or accidental relationship) 
between the coefficients of the XXI and XYZ terms that is responsible for the existence of this DFS. 



9 Note that our indices here differ somewhat from the notations in Sec. [VB, because there we considered either one- 
dimensional irreps, or mostly avoided explicitly indicating superpositions between copies of a given irrep. 
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